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This  is  a  compilation  of  various  parametric  conditions  for  the  critical 
ionization  velocity  (CIV)  process  to  occur  in  the  ionosphere. 
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Abstract — This  is  a  compilation  of  various  parametric  condi¬ 
tions  for  the  critical  ionization  velocity  (CIV)  process  to  occur  in 
the  ionosphere. 

I.  Introduction 

Tests  of  critical  ionization  velocity  (CIV)  in  the  labo¬ 
ratory  have  repeatedly  demonstrated  its  occurrence  with 
ease,  while  the  space  experiments  have  yielded,  at  best,  mixed 
results  (see,  for  example  [1]).  Whether  or  not  CIV  occurs  in 
space  and  how  easily  it  is  initiated  have  important  implications 
in  the  Space  Shuttle  environment  (see.  for  example,  [2]).  In 
order  to  demonstrate  CIV  clearly  in  space,  it  is  necessary  to 
consider  the  optima!  conditions  under  which  CIV  may  occur. 

Although  Alfven  originally  suggested  in  the  CIV  concept 
that  rapid  ionization  would  occur  when  the  critical  ionization 
velocity,  V'»,  is  exceeded,  transferring  this  concept  into  an 
actual  positive  space  experiment  requires  that  careful  attention 
be  paid  to  determining  the  details  of  such  experiments.  In  view 
of  the  increasing  number  of  CIV  space  experiments  conducted 
or  planned  in  recent  years,  it  is  useful  to  compile  the  main 
inequality  conditions  which  have  to  be  satisfied  for  CIV  to  be 
triggered.  This  paper  is  not  meant  to  be  an  exhaustive  survey, 
but  it  aims  at  discussing  the  most  important  and  salient  aspects 
of  CIV. 

II.  Beam  Velocity 

Alfven  (3],  [4]  postulated  that  when  the  relative  velocity. 
V,  between  a  neutral  gas  and  a  magnetized  plasma  exceeds  a 
critical  value,  V*,  rapid  ionization  occurs.  Alfven  [3].  [4]  gave 
a  formula  for  the  critical  ionization  velocity,  \\,  as 

V*  =  yj2ecplM.  (1) 

where  eo  is  the  ionization  energy  and  M  the  mass  of  a  neutral 
atom  or  molecule.  Subsequent  laboratory  CIV  experiments 
(e.g.  [5])  have  shown  that  rapid  ionization  occurs  at  values  of¬ 
ten  slightly  higher  than  V*.  Foimisano  et  at.  [6]  introduced  an 
empirical  efficiency  factor,  rj,  into  Alfven’s  formula  (eq.  (1)) 
to  account  for  energy  loss  by  various  (sometimes  intangible) 
mechanisms: 


where  r?  is  smaller  than  unity.  Equation  (2)  replaces  (1)  if  t) 
is  significant.  For  critical  velocity  ionization  to  occur,  it  is 
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Fig.  1,  Parametric  domain  of  CIV  (adapted  from  [9]) 
necessary  for  the  velocity  V  to  exceed  V*: 

V  >  v;.  (3) 

Papadopouios  [7]  pointed  out  that  if  the  velocity  T  is  so 
fast  that  it  exceeds  a  critical  value  1',.  electromagnetic  modes 
may  be  excited; 

i;  =  {1  - (4) 

where  3^  ~  nkT/B^  and  \'a  is  the  Alfven  velocity.  These 
electromagnetic  modes  tend  to  drain  away  energy  [8].  Indeed. 
Machida  and  Goertz  [9]  and  Goertz  ei  al.  [  lOj  have  obtained 
numerical  simulation  results  demonstrating  the  suppression  of 
CIV  when  V  >  V,  with  the  modified  two-stream  instability 
and  electron  heating  becoming  inefficient.  Fig.  1  shows  a 
parametric  domain  of  CIV  with  boundaries  of  V  >  L*  and 
V'  <  Vi,.  Laboratory  measurements  [11]  on  CIV  of  H;0 
have  confirmed  the  electromagnetic  effect  in  quenching  CIV 
(Fig.  2). 

Thus,  for  CIV  to  occur,  the  beam  velocity  criterion  is 

V;  >  V  >  W.  (5) 

III.  Townsend  s  Criterion 

Townsend’s  criterion  [12]  for  discharges  states  that,  in  order 
to  sustain  a  discharge,  an  electron  has  to  generate  at  least  one 
electron  during  its  residence  time  in  the  interaction  region.  In 
typical  CIV  space  experiments,  a  beam  or  cloud  of  neutral 
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Fig.  2,  Ionization  density  n  as  a  function  of  magnetic  field  B  in  a  latmratory 
CIV  experiment  (from  [Ilj). 

gas  is  injected  across  the  ambient  magnetic  field.  The  beam, 
or  cloud,  diverges  as  it  propagates.  The  ionization  time,  r,,  is 
given  by 

r,(x)  =  l/;n(j-)(T(iv)i>’,  (6) 

where  x  is  the  distance  from  the  exit  point,  rt  the  neutral 
gas  density,  a  the  cross  section  of  electron  impact  ionization, 
and  Vg  the  electron  velocity.  To  satisfy  Townsend's  criterion, 
the  ionization  time,  r,.  has  to  be  shorter  than  the  electron 
residence  time,  Tg.  In  the  literature,  two  residence  times  have 
been  considered.  TTiey  are  the  electron  transit  time  [2]  and 
the  cloud  transit  time  [13],  [14|.  These  two  versions  will  be 
discussed  as  Townsend’s  conditions  f  and  II  respectively  in 
the  following. 

Townsend’s  Condition  I 

Consider  a  conical  beam  of  neutral  gas.  The  neutral  density, 
n(x),  drops  as  the  beam  expands: 

n(i)  = -Vq/ [Wx*  tan"((i)/2)] .  (7) 

where  Nq  is  the  release  rate,  V  the  beam  velocity,  and  o  the 
cone  angle.  The  electron  transit  time,  Tc(x),  increases  with  x 
because  of  the  beam  divergence: 

Tgix)  —  2xtan{0/2)v~^  (8) 

Townsend’s  condition  in  this  case  is 

tJx)  >  Tiix).  (9) 

which  yields  a  critical  distance,  xj,  beyond  which  Townsend's 
condition  is  not  satisfied.  This  implies  that  CIV  can  only  occur 
within  a  distance  x  <  xi,  which  is  given  by 

xi  =  2iVocr/[V'7rtan(d/2)|.  (10) 

Fig.  3  shows  the  critical  distance  as  a  function  of  No-  Equation 
(6)  is  valid  for  a  steady-state  continuous  beam  only  or  for  a 
gas  cloud  of  length  d  as  long  as  d  >  xj. 

Consider  next  a  gas  cloud  of  length  d  with  d  <  xj.  The 
density,  njjx),  of  the  cloud  can  be  modeled  as 

nrf(x)  =  3A7[47r{Vr<)^] .  (11) 

where  N  is  the  total  number  of  neutrals  at  f  =  0,  V).  is 
the  radial  expansion  velocity,  the  cloud  length  d  =  Ij-f,  and 
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Fig.  3.  Townsend  distance  t\  in  conical  beam  mixiel  (from  j2j) 


the  distance  from  the  exit  point  is  x  =  \'t.  In  this  ca.se. 
Townsend’s  condition  (eq.  (9))  yields  a  critical  distance  xj: 

xi  =  Xd  =  (V/\;)(3.Vn:2T('  ^  (12) 

Fig.  4  shows  the  Townsend  critical  distance  computed  for 
CRIT  II  and  CRRES-Iike  experimental  situations. 

We  remark  that  both  of  the  Townsend  I  results  (10)  and 
(12),  are  independent  of  the  electron  velocity,  v,..  Therefore, 
even  if  the  electrons  are  subsequently  slowed  down  by  ion 
space  charge,  the  result  is  the  same.  However,  if  the  effective 
distance  traveled  by  the  electrons  becomes  longer,  due  lo 
scattering  for  example,  then  the  requirements  (10)  and  (12) 
would  become  less  stringent. 

Townsend's  Condition  II 

In  CIV,  electron  energization,  or  heating,  occurs  in  a  time 
interval  Tf{ .  As  a  neutral  cloud  propagates  forw  ard  across  the 
magnetic  field,  the  electrons  being  energized  are  more  or  less 
tied  down  by  the  field  lines.  Therefore  the  electrons  lag  behind 
the  ions,  which  arc  almost  unmagnetized  for  time  shorter  than  a 
quarter  of  the  ion  cyclotron  period.  If  the  cloud  dimension  d  is 
short,  the  electrons  may  lag  behind  the  cloud  with  insufficient 
heating.  Therefore,  the  following  condition  jl8]  is  necessary 
for  CIV: 


^'Ch‘  >  *• 


where 


=  ,hV 


(14) 
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Fig.  4.  Calculated  critical  distance  j-j  for  CRIT-II  and  CRRES  like  releases. 
The  one  for  CRIT-II  used  a  random  velocity  of  ±15  km/s  and  a  mam 
velocity  of  T  km/s.  The  one  for  CRRES  used  a  main  velocity  of  11  km/s 
and  a  random  velocity  of  ±1  km/s  (from  |16j). 


Fig.  5.  Pulse  expansion  in  the  Star  of  Lima  experiment.  The  points  are  the 
locations  of  pulse  heights  estimated  from  Torben's  data  [17J  (from  (18J). 

is  the  transit  time  of  the  cloud  across  a  magnetic  field  line. 
Fig.  5  illustrates  the  transit  of  an  expanding  cloud  in  a  CIV 
space  experimental  situation. 

For  example,  if  an  initially  (^-function  pulse  cloud  is  injected 
and  the  cloud  expands  according  to  rf  =  \\-t,  then  (13) 
requires 

t  >  (y/V^}rH  (15) 

or 

X  >  T2  =  (V'^/Vr)r//,  (16) 

which  implies  CIV  cannot  occur  at  distances  too  near  the  exit 
point.  Combining  Townsend  s  conditions  I  and  II,  we  find  that 
CIV  has  to  occur  between  the  critical  distances  xj  and  X2: 

(17) 


TABlfc  1 
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Before  closing  this  section,  two  remarks  will  be  gj\en 
The  first  remark  concerns  the  heating  time,  th  In  the  above 
discussion  (cqs.  {13)-(17)).  we  have  not  specified  ~n-  which 
depends  on  the  healing  situation.  If  the  electrons  are  heated 
from  an  initial  state  in  which  there  are  a  stationary  ambient 
plasma  and  .some  ion  seeding  mechanism  which  builds  up  an 
ion  beam  density  from  zero  as  considered  by  Brenning  jl4j 
and  Anxiis  (15],  is  extremely  sensitive  to  the  seeding  rate 
and  the  initial  electron  temperature.  We  denote  this  type  of 
Tfj  as  7^1.  If,  on  the  other  hand,  the  electrons  are  heated 
when  an  ion  beam  is  already  formed,  as  considered  by  Tanaka 
and  Papadopoulos  (31).  a  good  estimate  of  "h  is  IO/wth 
according  to  computer  simulations  (31).  W'e  denote  this  type 
of  Th  as  THi-  For  a  continuous  injection  as  considered  in  this 
paper,  it  is  likely  that  the  latter  situation  (3!)  is  more  relevant. 
In  the  rest  of  the  paper,  is 

The  second  remark  concerns  another  version  of  Townsend  s 
criterion.  Haerendel  [13]  used  ionization  time.  ",  in  place  of 
heating  lime,  in  (13).  This  places  much  milder  restrictions 
on  the  CIV  process,  and  will  not  be  considered.  For  experi¬ 
mental  designs,  it  is  more  cautious  to  u.se  a  restnciise  criterion 
(eq.  (13)).  Table  I  summarizes  the  combinations  discussed  in 
the  two  remarks. 

IV.  Magnetic  Fieut  Strength 

Brenning  [19]  proposed  on  theoretical  grounds  an  upper 
limit  to  8  (eq.  (18))  for  CIV  to  cKCur.  Equation  (IS),  to  utir 
knowledge,  has  not  yet  been  tested  experimentally: 

a,',  (n,.  i/dvr(  Z?i  >  I*  ■  (IS) 

Here  w.  and  .vv  are  the  plasma  and  electron  gyrofrequencies 
respectively,  rn  the  electron  mass.  M  the  ion  masses,  and  n, 
the  plasma  density. 

Brenning  [8]  pointed  out  that  (4)  gives  a  lower  limit 
of  magnetic  field.  Brenning  [19]  showed  experimentally  an 
absence  of  CIV  for  a  weak  magnetic  field.  The  absence  of 
CIV  in  Breniiing's  experiment  [19]  agrees  completely  with 
the  lower  B  field  limit  given  by  (4). 

Suppose  for  a  given  set  of  experimental  parameters,  the 
value  of  the  magnetic  field  satisfying  (18)  is  Di  and  that 
satisfying  (4)  is  Bn.  Then,  the  range  of  magnetic  field  B  for 
CIV  to  occur  is  limited  by 

Bx  >  B  >  «  (Un 

where  Bi  is  given  by  (18)  and  Bn  by  (4). 


Xl  >  T  >  X2- 


LAI  .\.ND  MURAD:  INEOUALITV  CONDITIONS 


n 


I 


ni 


-L  L  Z 


Fig.  6.  The  growth  rate.  !m  decreasing  as  kL  decreases,  in  a  simple 
density  profile  beam  mode!  (from  (22]). 


Fig.  t  CIS'  parametric  domain,  il  shows  that  WO  '  mas  noi  he  oplimai 
(from  (2|) 


V.  Beam  Width 

Large  beam  width  L  may  favor  CIV  m  two  ways.  Firstly, 
larger  beam  width  implies  that  the  residence  time  of  an 
electron  is  longer; 

r,  =  L/r,.  (20) 

where  is  the  electron  velocity  along  the  magnetic  field 
direction.  Longer  in  turn,  implies  higher  probability  of 
electron  impact  ionization  in  the  beam  and  hence  better  chance 
for  Townsend’s  condition  I  (eq.  (9))  to  be  satisfied.  Besides, 
escape  of  hot  electrons  implies  not  only  a  reduction  in  the 
number  of  electrons  but  also  energy  loss.  Therefore,  the  wider 
the  beam,  the  better  the  chance  for  CIV.  Secondly,  Newell 
[20]  and  Kelley  et  al.  [21]  have  suggested  that  in  order  for  the 
modified  two  stream  instability  to  be  efficient,  it  is  necessary 
for  the  beam  to  be  wide  enough  to  accommodate  the  parallel 
component  of  the  lower  hybrid  wavelength: 

L>\.  (21) 

A  simple  theoretical  square  density  model  [22]  has  been  found 
to  yield  a  monotonically  decreasing  function  of  the  instability- 
growth  rate  Im  w/w,,  where  a.-,  is  the  ion  plasma  frequency,  as 
a  function  of  decreasing  beam  width  L  (Fig.  6).  The  theoretical 
result  agrees  with  the  suggestions  of  Newell  [20]  and  Kellev 
et  al.  [21]. 

In  typical  CIV  space  experiments,  the  value  of  A||  is  of  the 
order  of  km  [21].  For  a  diverging  beam  with  a  longitudinal 
velocity  of  8  km/s  and  a  radial  velocity  of  1  km/s.  for 
example,  the  efficiency  for  CIV  is  reduced  by  this  effect  in 
the  early  1  second  period  or  8  km  distance.  Conceivably,  the 
proposed  CIV  process  in  the  primordial  solar  system  [3]  can 
satisfy  (21)  because  the  dimension  L  was  astronomically  large. 

We  remark  that  violation  of  the  condition  L  >  A||  (eq.  (21)) 
does  not  preclude  CIV  although  it  would  reduce  the  instability 
growth  rate. 

Finally,  we  remark  that  there  is  no  known  upper  limit  on 
the  beam  width  L  for  CIV  to  occur.  The  efficiency  of  CIV 
increases  monotonically  with  the  width  of  the  beam.  If  CIV 
can  occur  in  the  spacecraft  environment,  the  environment  of 


the  space  station  with  its  large  size  (order  of  km)  may  be  a 
more  likely  place  for  CIV  to  occur. 

VI.  Pitch  A.ngle 

We  discuss  here  two  aspects;  1)  electron  escape  and 
2)  electron  energization.  In  the  first,  we  consider  the  increase 
in  the  effective  distance.  L,  for  an  electron  traveling  along  B. 
or  in  the  plane  of  B  and  E  x  B  The  magnetic  pitch  angle  H 
of  a  beam  need  not  be  90°: 

Ue}  =  L/^inft.  (22) 

where  L  =  L{9  =  90°  ).  However,  the  beam  kinetic  energy 
E  available  for  energizing  the  instability  waves  and  electrons 
would  be  reduced  accordingly: 

=  90°)  sin"  (23) 

where  E{0  =  9O°}  =  {l/2)A/r-.  Since  -  (eq.  (20))  has  to 
be  greater  than  r,  (eq.  (9)).  and  E{9)  has  to  be  greater  than 
eo  (eq.  (1)).  one  obtains  from  (1),  (9),  (20).  (22).  and  (23)  an 
inequality  [2]: 

Lan  >  siri^  >  V,/T’.  (24) 

where  L,  a.  n.  V'*,  and  1'  have  been  defined  already.  Fig.  7 

shows  a  parametric  domain  of  CIV  as  limited  by  this  aspect. 

With  regard  to  electron  energization.  Goertz  etal.  [23]  have 
emphasized  that,  in  the  heating  process  of  electrons  by  means 
of  the  modified  two-stream  instability,  B  -  90°  may  not 
always  be  optimum  for  two  reasons.  The  first  is  that  resonance 
w'ave- particle  heating  at  an  electron  tail  can  only  heat  a  small 
electron  population  to  even  higher  energy.  It  is  better  to  heat 
just  below  the  tail,  where  the  population  is  higher,  so  that 
more  electrons  will  be  heated.  The  electrons  heated  in  such  a 
way  then  become  hot  electrons  as  they  migrate  into  the  tail 
distribution.  Therefore,  for  a  given  beam  energy,  the  pitch 
angle  0  should  be  adjusted  so  that  the  wave  phase  velocity  is 
just  below  the  tail  in  order  to  achieve  optimum  heating  [24], 
The  second  reason  is  that  the  cross  section  o  of  ionization 
features  a  maximum  around  70  to  100  eV  for  most  elements. 
If  the  beam  has  too  much  kinetic  energy,  the  electrons  may  be 
heated  to  an  energy  too  high  for  optimum  ionization.  Adjusting 
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Ionization  Fraction  Ba 

Fig.  8  Crilical  ionization  velocity  in  a  mixture  of  Sr  and  Ba  (from  (25]) 

6  may  yield  a  hot  tail  at  the  peak  a  range.  However,  we  must 
remark  that  {<tv^)  is  often  approximately  constant  in  many 
cases.  Thus,  this  reason  is  less  important.  Also,  in  order  to 
achieve  maximum  ionization  rate  for  a  given  available  energy 
Vr,  one  maximizes  {av^)/W,  where  W  =  {l/2)Tnr;  for 
monoenergetic  electrons.  The  maximization  condition  gives 
dajdv  —  a/i',  which  requires  energies  below  the  maximum 
of  (7(f). 

VII.  Mixture  of  Species 

Anxas  (5]  reported  experimental  results  which  show  only 
one  critical  ionization  velocity  in  a  binary  mixture  of  species, 
while  Lai  et  al.  (25]  reported  a  theoretical  analysis  of  a  general 
mixture  of  species.  In  a  mixture  of  species  with  a  given 
velocity  distribution  f{V)  and  various  masses  and  ionization 
potentials,  some  species  may  act  as  energy  sources  while 
others  as  energy  sinks.  Depending  on  the  parameter  domains, 
it  is  possible  to  have  single  or  multiple  species  undergoing 
CIV.  A  generalized  Alfven  CIV  formula,  which  represents  the 
energy  budget,  has  been  suggested  [25]: 

J dV[{l/2)M,V^  -  e<p.]/(V-)  >  0,  (25) 

where  i  labels  the  species,  V  the  neutral  velocity,  and  v  the 
fractional  ionization.  Fig.  8  shows  this  energy  budget  for  a 
binary  mixture  of  Sr  and  Ba. 

Depending  on  the  velocity  V  and  the  extent  of  the  hot  tail  of 
electrons,  single  and  multiple  critical  ionization  velocities  may 
exist  [25].  When  the  number  of  species  is  1,  (25)  reduces  to 
Alfven 's  formula  (eq.  (1));  when  the  number  is  2,  (25)  reduces 
to  Axnas’s  formula  [5]. 

To  illustrate  this  treatment,  we  consider  the  Space  Shuttle 
exhaust,  which  is  a  mixture  of  several  species  (2].  The 
main  species  are  H2O,  Ni,  H2,  CO,  and  CO2,  their  mole 
fractions  being  0.33,  0.31,  0.16,  0.13,  and  0.042  respectively 
(25}.  In  the  simplest  case,  all  species  travel  at  the  same 
velocity  V.  The  critical  ionization  velocity  calculated  using 
(25)  [25]  is  very  near  the  classical  Alfven  critical  velocity  of 
H20(V'*  «  il.6  km/s).  As  another  example,  a  square  shaped 
velocity  distnbution  f{V)  =  11.2  ±  AV',  would  require 


NEON  BEAM  velocity^!  !  v, 
INITIAL  NEUTRAL  DENSITY: 


Fig.  9.  Effect  of  radiative  lifetime  of  the  met3.vlable  state  on  the  onset  of 
ionization  (from  (27)) 

AV  >  ±5.2  km/s  in  order  to  satisfy  the  energy  budget  and 
achieve  CIV. 

VIII.  Metastable  States 

As  discussed  earlier,  hot  electrons,  which  are  in  the  tail 
of  the  Maxwellian  distribution,  are  much  less  abundant  than 
warm  electrons,  which  are  just  below  the  tail.  When  the  energy 
budget  is  tight,  these  electrons  alone  mav  have  difficulty 
in  triggering  and  sustaining  a  CIV  process.  If  the  lower 
energy  (and  more  abundant)  electrons  can  also  participate  in 
ionization,  it  would  be  easier  to  achieve  CIV.  In  discharge  tube 
studies,  the  importance  of  metastable  states  in  ionization  and 
gas  breakdown  has  been  established  (e.g.  [26]).  In  a  two-step 
ionization  process  via  the  metastable  state,  each  step  requires 
an  energy  less  than  eo,  the  ground-state  ionization  energy. 
Therefore,  metastable  states  may  foster  CIV,  especially  under 
tight  budget  situations.  Lai  et  al.  [27]  and  McNeil  et  al.  [28] 
have  demonstrated  that  metastable  states  can  indeed  help  CIV 
provided  that  the  metastable  state  lifetime.  r\;.  is  long  enough. 
For  example,  in  a  simple  model  of  neonlike  CIV  discharge 
[27],  inclusion  of  metastable  states  would  make  at  least  SO'T- 
faster  ionization  if  r.v/  >Sxl0~^s  (Fig.  9). 

To  estimate  the  critical  r\/.  one  equates  the  metastable  state 
lifetime,  r.w,  to  the  ionization  time  per  atom  for  the  transition 
(m  ^  i)  from  a  metastable  state  to  the  ionizing  state.  One 
obtains 

Tsr  >  (26) 

where  ric  and  are  the  electron  density  and  velocity,  and 
<7rni  the  cross  section  of  the  transition  (w  —  7).  Larger  values 
of  lifetime  t^{.  electron  velocity,  electron  density,  and  cross 
section  17^1  would  favor  metastable  state  ionization. 

IX.  Mass  Loading 

Haerendel  [13],  [29]  pointed  out  that  the  beam  ions,  travel¬ 
ing  with  a  velocity  V  .  inevitably  impart  momentum  to  the 
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Fig.  10.  A  conceptual  picture  of  momentum  coupling  of  CIV  in  space.  The 
beam  travels  for  a  distance  R  while  the  Alfven  waves  a  distance  D 


ambient  ions.  As  a  result,  the  ambient  ions,  together  with 
the  corresponding  magnetic  flux  tube,  love  forward,  reducing 
the  relative  velocity,  14,  between  the  beam  and  the  ambients. 
The  phenomenon  was  called  mass  loading  [13],  [29].  A  simple 
estimate  of  the  critical  beam  density,  above  which  mass 
loading  becomes  significant,  can  be  obtained  as  follows. 

Consider  a  beam  of  dimension  R.  Following  Haerendel  [13], 
[29],  we  rssume  that  the  length  D  of  the  magnetic  flux  tube 
involved  in  time  T  is  given  by 

D  =  TV^.  (27) 

where  is  the  Alfven  velocity,  and  T  the  time  required  for 
tlie  beam  to  traverse  the  flux  tube  (Fig.  10)  such  that 

R  =  TV.  (28) 

If  the  flux  tube  travels  with  a  fraction  Sa  of  the  beam  velocity 
V ,  we  equate  tne  momenta  lost  by  the  beam  and  gained  by 
the  flux  tube: 


MbnbR{\  -  Sb)  =  2Man^R(V\/V)sa.  (29) 


where  the  subscripts  a  and  b  denote  ambient  and  ion  beam 
respectively,  Sb  is  the  fractional  slowdown  of  the  initial  beam 
velocity,  and  the  2  is  for  two  flux  tubes  (upper  and  lower). 
From  (29),  we  obtain  a  critical  beam  density  rib'. 


Tiff 


=  2na 


Mb 


V 


Sg 

1-56 


(30) 


If  the  flux  tube  and  the  injected  cloud  gets  a  common  velocity 
which  is  a  fraction  s  of  the  initial  beam  velocity,  s  =  Sg  —  Sb- 
We  obtain  a  critical  beam  density  at  which  s  =  Sg  =  Sb  =  0.5 
(i.e.,  the  injected  beam  is  braked  to  half  its  initial  velocity): 


or 


Tib  —  271^ 


MgVA 

MbV 


B  I  rig  Mg 

~  MbV  V  7T 


(31) 


(32) 


approximately  constant  while  the  ambient  plasma  density. 
varies  in  day  and  night  or  solar  cycles.  (30)  or  (32)  imply 

Tib  >.  yN,,  . 

which  means  higher  allows  higher  critical  beam  density  n* 

X.  AN  E.k,\.mple.  of  Multiple  Ineoualities 

As  an  exercise,  we  consider  «  composite  case  in  which 
several  of  the  instabilities  discussed  are  involved.  We  consider 
the  scenario  of  a  hypothetical  space  CIV  experiment,  in  which 
a  xenon  gas  is  released  from  a  satellite  traveling  at  7.4  km  s 
across  the  ambient  magnetic  field  in  the  ionosphere.  We 
assume  the  following  initial  conditions: 

Gas  pressure,  p,  in  storage  chamber  =  10  atmospheres; 

Gas  temperature.  T,  in  storage  chamber  -  2T.3'  K, 

Radius  of  exit  aperture  =  0.25  cm: 

Beam  half  cone  angle  0/2  =  30'; 

Gas  is  pulsed  in  a  sequence  of  10  s.  1  s.  and  a  shorter 
one. 

The  first  inequality  to  be  used  is  (5).  At  about  5(X)  km  alti¬ 
tude,  the  Alfven  velocity,  1 4.  for  oxygen  is  about  7{X)  km/s. 
and  i3e  Js;  2  X  10“^  (30).  Thus.  (4)  gives  V  ^  I'a  % 
700  km/s  for  oxygen  or  about  240  km/s  for  Xe.  As  long  as 
the  xenon  beam  velocity  is  between  240  km/s  and  4.225  km/s 
approximately,  we  expect  the  velocity  requirement  (eq.  (5)) 
for  CIV  to  be  fulfilled. 

We  now  use  the  second  inequality  (eq.  (17))  From  gas 
kinetic  considerations,  the  initial  gas  velocity.  V.  is  given  in 
the  spacecraft  frame  of  reference  by 

V  =  (8ifcT'/-.\/)’  "  =  210  m/s.  (.33) 

where  Vf  is  the  molecular  weight.  The  initial  mass  density  p 
is  a  function  of  pressure  p; 

p  =  pr-/3.  (34) 

The  mass  release  rate  R  is  given  by  /?  =  I  '.4p.  w  here  .4  is  the 
aperture  area.  The  molecular  release  rate  .V„  =  R/M.  With 
these  conditions,  the  first  Townsend  distance  xi  (eq.  (10))  in 
the  spacecraft  frame  of  reference  becomes 

X]  =  I'-iiAp/kTicr.  (35) 

Remarkably.  (35)  is  independent  of  the  molecular  weight  M 
of  the  gas  species.  In  this  sense,  (38)  is  a  universal  formula 
for  the  “Townsend  I”  distance  X}.  Taking  a  cross  section 
a  _  10“*®  cnr.  (38)  gives  xi  =  66  m  from  the  exit  point. 
For  example,  if  there  are  several  gas  releases  with  different 
molecular  masses,  the  universal  formula  (eq.  (35))  yields  the 
same  “Townsend  I”  distance  for  all  releases  provided  that  the 
cross  section  a  is  the  same. 

For  a  pulse  of  1  second  duration,  the  pulse  length  rl  is 
0.21  km  from  the  exit  point.  Assuming  that  the  ambient 
magnetic  field  is  stationary  in  the  space  frame,  Ihe  time 
for  a  field  line  to  sweep  through  the  pulse  is  given  by 

T,/  =  0.21 /( 7.4  +  0.21)  ==  27  ms. 


Remarkably,  mass  loading  implies  that  higher  ambient  plasma  Suppose  the  electron  heating  time,  r//,  is  given  by  .30.  x/  // 
density  favors  CIV.  If  the  altitude  and  therefore  M„  remain  [31}  and  x/,H(  Ba)  is  1034  Hz  [30];  we  have  r/^i  Xo)  s: 
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TABLE  H 

ISEOUAUTIES 


Paiameler 

Inequality 

Equati'’” 

Beam  velocity,  I ' 

t ,  >  t  >  u. 

Distance,  x 

Xt  >  X  >  X  '.> 

p 

Magnetic  field.  B 

Bi  >  B  >  0j 

IV 

Beam  width.  L 

L  >  .V 

21 

Beam  angle.  0 

Ltrn  >  %iu0  >  1  .,'T' 

:4 

Fractional  ionization,  i', 

53,  If,  /dU[(l/2)T/,U-  -^o,i/(r  :  >  U 

25 

Metastable  lifetime,  r.\/ 

7.-.J  >  ) 

:6 

Beam  ion  density,  rtf, 

ns  <  2r,  jM^/XhlO  A /I'! 

31 

27  ms,  which  happens  to  equal  rj.  Thus,  the  1  second  pulse 
may  marginally  satisfy  the  Townsend  II  criterion  (eq.  (13)). 
The  shorter  pulse  ejected  later  will  not  satisfy  (13)  and 
therefore  will  not  undergo  CT'. 

Consider  next  the  plasma  density  «  2.5  x  lO^/cc;  here 
(18)  yields  Bi  =5  800  G  for  Xe*.  In  the  ionosphere,  B  is 
typically  0.22  G.  Therefore  there  is  no  difficulty  in  satisfying 
Bi  >  B.  As  for  the  second  inequality  of  (4),  we  observe  that 
since  0  ^2x  10“^  in  the  ionosphere  [30],  it  takes  an  ambient 
magnetic  field  five  orders  of  magnitude  weaker  to  have  a 
significant  effect.  Therefore,  the  magnetic  field  inequality  (19) 
is  satisfied  easily  in  CIV  space  experiments. 

In  order  to  satisfy  (21),  the  conical  pulse  has  to  travel  for 
at  least  6  s  before  its  transverse  dimension  L  reaches  1  km 
approximately.  Prior  to  6  s,  the  dimension  L  is  too  small 
for  the  full  efficiency  of  the  modified  two-stream  instability 
growth  rate  to  take  place.  Failure  to  satisfy  (21)  does  not 
necessarily  rule  out  CIV  although  the  efficiency  is  reduced. 

To  illustrate  the  utility  of  (24),  suppose  we  wish  to  have  CIV 
at  X  =  70  m,  which  is  located  just  beyond  the  Townsend  I 
distance  (ii  ss  66  m),  which  is  obtained  by  assuming  perpen¬ 
dicular  propagation  (i.e.,  6  —  90°).  Equation  (7)  gives  n(x)  = 
1.2  X  10^^/cc.  The  transverse  distance  at  i  is  L  =  80.8  m.  The 
first  inequality  of  (24)  gives  Lcrn  ~  0.972  >  sin  restricting 
6  to  be  below  76°,  the  value  for  cr  being  taken  as  10“^®  cm^. 
The  second  inequality  of  (24)  requires  that  the  beam  velocity 
V  exceeds  4.35  km/s  to  achieve  CIV  for  Xe.  This  value  is 
slightly  higher  than  the  normal  xenon  CIV  value  (4.23  km/s). 

Returning  to  the  scenario,  we  now  consider  metastable 
states.  Xenon  has  two  metastable  states,  one  of  which 
(6s[3/2]2)  has  a  lifetime  tm  of  150  s  at  8.32  eV  [32]. 
Ionization  from  a  metastable  state  may  help  CIV.  However, 
if  a  metastable  state  is  not  ionized  fast  enough,  it  may  decay 
to  a  lower  state  with  loss  of  energy  in  the  form  of  radiation. 
To  estimate  a  critical  electron  density  rie,  we  use  (26),  which 
gives 

Tie  >  ll(TMarmV^).  (36) 

Using  the  data  of  Cartwright  et  al.  [33],  the  cross  section  Omi 
has  been  modeled  [27]  as  about  10~^°  cm^  in  the  electron 
energy  range  =  o  io  25  eV.  Equation  (36)  requires  an 
electron  density  >  4  x  10“*  or  2.3  x  10'*/cc  for  £,  =  5 
or  25  eV,  respectively,  for  the  metastable  state  to  play  a 
significant  role  in  CIV  of  xenon.  Such  electron  densities  can  be 
easily  met  in  the  ionosphere,  even  at  an  altitude  near  800  km. 


As  an  estimate  of  mass  loading,  we  take  n„  s:  10\.  <  c.  ^ 
700  km/s  [34],  Ma{0)/Mf,(Xe)  =:  1/9,  and  V  %  7,8  km/s. 
Taking  Sg  =0.3  and  .sj,  =  O.I.  so  that  1)  =  (1  -  -  .</,  ;.r  = 

0.6V  =  4.68  km/s  in  (29).  one  obtains  rii,  a  6.6  x  10"/ct. 
Since  the  value  of  V,  for  xenon  is  about  4.2  km/s.  the  value 
of  1/  is  close  to  V,.  In  view  of  the  possible  presence  of 
other  intangible  factors  of  energy  loss,  it  is  possible  that  the 
efficiency  of  CIV  is  already  near  zero  at  )/•  Therefore,  we 
expect  an  upper  bound  of  the  ionization  level  of  the  xenon 
CIV  gas  release  to  be  about  6.6  x  10°/cc. 

XI.  Conclusion  and  Discussion 

The  purpose  of  this  paper  has  been  to  compile  some 
inequality  conditions  for  CIV  to  occur,  which  are  summarized 
in  Table  II,  and  to  explore  how  these  inequalities  limit  the  CIV 
conditions.  For  CIV  to  occur,  the  beam  velocity  I  ’  must  exceed 
a  critical  value,  V,,  but  too  high  a  velocity  V(>  I/,)  may 
lower  the  efficiency  (eq,  (5)).  For  a  diverging  neutral  beam 
injected  into  space,  Townsend’s  criteria  1  and  I!  (eq.  (17)) 
impose  two  critical  distances,  X]  and  X2;  CIV  can  occur  only 
between  these  critical  distances.  The  ambient  magnetic  field  in 
the  ionosphere  is  well  within  the  inequality  bounds  (eq.  (19)). 
For  maximum  growth  of  lower  hybrid  modes,  the  parallel 
wavelength  Aij  has  to  fit  into  the  dimension  of  the  interaction 
region  (eq.  (21)).  This  dimension  requirement,  together  with 
Townsend’s  criterion  II,  implies  that  CIV  is  likely  ruled  out  ai 
distances  (of  the  order  of  meters  typically)  too  near  the  beam 
origin.  Beam  angle  &  =  90°  may  not  be  optimal  (eq.  (24));  a 
smaller  angle  may  allow  a  longer  parallel  electron  trajectory 
in  the  interaction  region.  For  a  mixture  of  species,  the  critical 
ionization  velocity  is  not  given  by  the  lowest  valued  species 
but  by  an  energy  budget  equation  (eq.  (25)).  Metastable  states 
can  foster  CIV  if  the  lifetime  of  the  states  is  long  and  the 
electron  density  is  high  (eq.  (29)).  When  the  beam  ionization 
is  too  high,  mass  loading  can  reduce  the  efficiency  of  CIV 
(eq.  (30)),  This  imposes  an  upper  limit  of  ionization  on  beams 
undergoing  CIV  in  the  ionosphere. 

Although  these  inequalities,  obtained  by  using  simple  mod¬ 
els,  may  be  too  stringent,  they  are  helpful  in  providing 
guidelines  for  space  experiments  designed  to  achieve  (or 
avoid)  CIV.  We  have  not  included  other  chemical  aspects 
such  as  molecular  dissociative  recombination,  ion-molecule 
reactions,  and  associative  ionization.  These  processes  may  in¬ 
fluence  the  progress  of  CIV,  but  do  not  affect  in  a  fundamental 
way  the  initiation  of  CIV, 
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